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Abstract 

Two-jet cross sections in deep inelastic scattering at HERA are calculated 
in next-to-leading order. The QCD corrections are implemented in a new 
ep — > n jets event generator, MEPJET, which allows to analyze arbitrary 
jet definition schemes and general cuts in terms of parton 4-momenta. First 
results are presented for the JADE, the cone and the kx schemes. For the 
T^-scheme, disagreement with previous results and large radiative corrections 
and recombination scheme dependencies are traced to a common origin. 
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Deep inelastic scattering (DIS) at HERA is a copious source of multi-jet events. Typical 
two-jet cross sections^ are in the 100 pb to few nb range and thus provide sufficiently high 
statistics for precision QCD tests Topics to be studied include the determination 

of Qisifiji) over a wide range of scales, measurement of the gluon structure function (via 
19 ^ QQ)-! ci-iid the study of internal jet structure. Clearly, next-to-leading order (NLO) 
QCD corrections are mandatory on the theoretical side. The dijet cross section, for example, 
is proportional to ^^(/iR) at leading order (LO), thus suggesting a direct measurement 
of the strong coupling constant. However, the LO calculation leaves the renormalization 
scale /i/j undetermined. The NLO corrections substantially reduce the renormalization and 
factorization scale dependencies which are present in the LO calculations and thus reliable 
cross section predictions in terms of asijnz) are made possible. NLO corrections in jet 
physics imply that a jet (in a given jet definition scheme) may consist of two partons. Thus 
first sensitivity to the internal jet structure is obtained, like dependence on the cone size 
or on recombination prescriptions. This will be of particular importance in the subsequent 
discussion. 

In this Letter, we investigate NLO QCD corrections to two-jet production in various 
jet algorithms at HERA energies. Previous calculations were limited to a JADE type 
algorithm, called l^-scheme below. In addition, approximations were made to the matrix 
elements which, as we will show, are not valid in large regions of phase space. Our calculation 
uses the s^m technique of Giele and Glover [^. This technique considerably simplifies the 
structure of NLO QCD corrections to hadronic processes and has already been applied to 
the calculation of NLO jet cross sections at LEP and the Tevatron @J^. The calculation 
is based on a QED-type factorization of soft gluon singularities and the use of universal 
"crossing functions" 0. This allows for a modular approach to NLO calculations. Our 
calculation of dijet production in DIS is implemented as a full NLO Monte Carlo program, 
MEPJET. The essential benefit of the Monte Carlo approach is that all hard phase space 
integrals are performed numerically. This makes the implementation of arbitrary cuts and 
jet algorithms a relatively easy task. 

Let us briefiy discuss the essential features of our calculation. Further details will be 
given in a subsequent publication. Deep inelastic electron proton scattering with several 
partons in the final state, 

e"(/) + p{P) e^(/') + proton remnant(p,.) -|- parton l(pi) + . . . -|- parton n{pn) iX) 

proceeds via the exchange of an intermediate vector boson \^ = 7, Z. In the following, Z- 
exchange will be neglected. We denote the momentum of the virtual photon, 7*, by g = / — 
its absolute square by Q^, the square of the final hadronic mass by W"^ = {P + q)'^ and use 
the standard scaling variables x = Q"^/ (2P ■ q) and y = P ■ q/ P ■ I. The general structure of 
the n-jet cross section in DIS is given by 

da'''"^[n - ]et] =J2 f dr] /a(r/,/i|) da^{p = rjP, asifil), fxl, fxl) (2) 



^In the following the jet due to the beam remnant is not included in the number of jets. 
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where the sum runs over incident partons a = q,q,g which carry a fraction t] of the pro- 
ton momentum, a"" denotes the partonic cross section from which colhnear initial state 
singularities have been factorized out at a scale fip and implicitly included in the scale de- 
pendent parton densities faiV: f^'p)- Born approximation, the subprocesses 7* + g 1 + 9, 
'J* + q q + g, and 7* + (? — ^ Q + q contribute to the two-jet cross section. At O^a^) the real 
emission corrections involve 'y* + q^q + g + g, Y + q^q + q + q^ l*+g^q + q + g and 
analogous anti-quark initiated processes. The corresponding cross sections are calculated 
by numerically evaluating the tree level helicity amplitudes as given in Ref. Similarly, 
the finite parts of the one-loop amplitudes are obtained by crossing the one-loop results for 
e^e^ qqg from Ref. 0. The numerical evaluation of helicity amplitudes has the advan- 
tage that the full spin structure is kept and, therefore, the MEPJET program allows for the 
calculation of all possible jet-jet and jet lepton correlations in NLO. In addition the NLO 
corrections for polarized electron on polarized proton scattering become available. 

In the one-loop amplitudes the ultraviolet divergencies are removed by MS renormal- 
ization which introduces a dependence on the renormalization scale /i/j. Infrared as well 
as collinear divergencies associated with the final state partons are cancelled against cor- 
responding divergencies of the one-loop contributions (see below). The remaining collinear 
initial state divergencies are factorized into the bare parton densities introducing a depen- 
dence on the factorization scale fip. In order to handle these singularities we follow Ref. 
and use the technique of universal "crossing functions". Starting from the results of the 
NLO calculation with all partons in the final state, i.e. e^e~ — n + 1 jets, where no such 
singularities occur after adding real and virtual contributions, the "crossing functions" con- 
tain the convolutions of the parton distribution functions with Altarelli-Parisi kernels and 
in addition take into account the crossing of a final state cluster to the intitial state, within 
the Smin cone as defined below. Note that all "plus" prescriptions are absorbed into the 
numerical evaluation of these crossing functions. 

The 3-parton final states need to be integrated over the entire phase space, including the 
unresolved regions, where only two jets are reconstructed according to a given jet definition 
scheme. In order to isolate the infrared as well as collinear divergencies associated with 
these unresolved regions the resolution parameter Smm is introduced. Soft and collinear 
approximations are used in the region where at least one pair of partons, including initial 
ones, has Sij = 2pi ■ pj < Smin and the soft and/or collinear final state parton is integrated 
over analytically. Adding this soft+collinear part to the virtual contributions gives a finite 
result for, effectively, 2-parton final states. In general this 2-parton contribution is negative 
and grows logarithmically in magnitude as s^m is decreased. This logarithmic growth is 
exactly cancelled by the increase in the 3 parton cross section, once s^m is small enough 
for the approximations to be valid. The integration over the 3-parton phase space with 
Sij > Smin is done by Monte-Carlo techniques (without using any approximations). Since, at 
each phase space point, the parton 4-momenta are available, the program is flexible enough 
to implement any jet definition algorithms or to impose arbitrary kinematical resolution 
and acceptance cuts on the final state particles. This is an essential advantage over existing 
programs such as DISJET |^. 

A powerful test of the numerical program is the Smin independence of the final result. 
Fig. 1 shows the inclusive dijet cross section as a function of Smin for the jet algorithms 
to be defined below. As mentioned before, Smin is an arbitrary theoretical parameter and 
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FIG. 1. Dependence of the inclusive two-jet cross section in the kT, cone, JADE, and the 
W^-scheme on Syyiiyi^ the two-parton resolution pciranieter. Partons are recombined in the -E'-schenie. 
Error bars represent statistical errors of the Monte Carlo program. For the fairly soft jet definition 
criteria described in the text, Smin independence is achieved for Smin < 0.1 GeV^ 

any measurable quantity should not depend on it. One observes that for values smaller 
than 0.1 GeV^ the results are indeed independent of Smin- The strong Smin dependence of 
the NLO cross sections for larger values shows that the soft and coUinear approximations 
used in the phase space region Sij < Smin are no longer valid, i.e. terms of 0{smin) and 
C(smmlnSm,m) becomc important. In general, one wants to choose Smin as large as possible 
to avoid large cancellations between the virtual+collinear+soft part {sij < Smin) and the 
hard part of the phase space (sjj > Smin)- Note that factor 10 cancellations occur between 
the effective 2-parton and 3-parton final states at the lowest Smin values in Fig. 1 and hence 
very high Monte Carlo statistics is required for these points, s^m independence is achieved 
at and below Smin = 0.1 GeV^ and we choose this value for our further studies. 

For these numerical studies, the standard set of parton distribution functions is MRS set 



D-' ||TT|. We employ the two loop formula for the strong coupling constant with = 230 



MeV, which is the value from the parton distribution functions. The value of is matched 
at the thresholds iir = rUq and the number of flavors is fixed to = 5 throughout, i.e. 
gluons are allowed to split into five fiavors of massless quarks. Unless stated otherwise, 
the renormalization scale and the factorization scale are set to — — 1/2 X^i Pxi^^ 
where denotes the magnitude of the transverse momentum of parton i in the Breit 
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TABLE I. Two-jet cross sections in DIS at HERA. Results are given at LO and NLO for the 
four jet definition schemes and acceptance cuts described in the text. 





two-jet 


two-jet exclusive 


two-jet inclusive 




LO 


NLO 


NLO 


cone 


1107 pb 


1047 pb 


1203 pb 


kT 


1067 pb 


946 pb 


1038 pb 


W 


1020 pb 


2061 pb 


2082 pb 


JADE 


1020 pb 


1473 pb 


1507 pb 



frame. A running QED fine structure constant <y{Q^) is used. The lepton and hadron 
beam energies are 27.5 and 820 GeV, respectively. A minimal set of kinematical cuts is 
imposed on the initial virtual photon and on the final state electron and jets. We require 
40 GeV^ < g2 < 2500 GeV^, 0.04 < ?/ < 1, an energy cut of E{e') > 10 GeV on the scattered 
electron, and a cut on the pseudo-rapidity r] = — lntan(^^/2) of the scattered lepton and jets 
of \rj\ < 3.5. In addition jets must have transverse momenta of at least 2 GeV in both the 
lab and the Breit frame. 

Within these general cuts four different jet definition schemes are considered for which 
we have chosen parameters such as to give similar LO cross sections (see Table I). In the 
ly-scheme the invariant mass squared Sij = {pi + Pj)^ is calculated for each pair of final 
state particles (including the proton remnant). If the pair with the smallest invariant mass 
squared is below i/cutW^, the pair is clustered according to a recombination scheme. Unless 
stated otherwise, and for all jet algorithms, we use the i?-scheme to recombine partons, i.e. 
the cluster momentum is taken as pi +Pj. This process is repeated until all invariant masses 
are above ycutW^. The resolution parameter i/cut is fixed to 0.02. 

The experimental analyses in are based on a variant of the IV-scheme, the "JADE" 
algorithm. It is obtained from the H^-scheme by replacing the invariant definition Sij = 
{Pi + Pj)^ by Mfj = 2EiEj(l — cos 9ij), where all quantities are defined in the laboratory 
frame. Neglecting the explicit mass terms pf and p| in the definition of M^^- causes substantial 
differences in jet cross sections between the W and the JADE scheme. 

In the cone algorithm (which is defined in the laboratory frame) the distance AR = 



{Ariy + (A0)^ between two partons decides whether they should be recombined to a single 
jet. Here the variables are the pseudo-rapidity t] and the azimuthal angle 0. We recombine 
partons with AR < 1. Furthermore, a cut on the jet transverse momenta of prij) > 5 GeV 
in the lab frame is imposed in addition to the 2 GeV Breit frame cut. 

For the algorithm (which is implemented in the Breit frame), we follow the description 
introduced in Ref. [jl2[. The hard scattering scale is fixed to 40 GeV^ and i/cut = 1 is 
the resolution parameter for resolving the macro-jets. In addition, jets are required to have 
a minimal transverse momentum of 5 GeV in the Breit frame. 

With these parameters, one obtains the two-jet cross sections of Table |. While the 
higher order corrections in the cone and kx schemes are small, very large corrections appear 
in the H^-scheme. In addition, the large effective ii'-factor (defined as i^' = (Jnlo/<^lo) 
of 2.04 (2.02) for the two-jet inclusive (exclusive) cross section in the VT-scheme depends 
strongly on the recombination scheme which is used in the clustering algorithm (see below). 
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FIG. 2. Dependence of a) the two-jet inclusive and b) the two-jet exclusive cross section in 
the cone scheme on the renormalization and factorization scale factor ^. The solid curves are for 

= /^F = ? Pt(0)^' while for the dashed curves only £,r = C varied but = 1/4 is fixed. 
Choosing the photon virtuality as the basic scale yields the dotted curves, which correspond to 
/^fl — l^^F ~ ^ Q"^- Results are shown for the LO (lower curves) and NLO calculations. 



A good measure of the improvement at NLO is provided by the residual scale dependence 
of the cross section. We have considered scales related to the scalar sum of the parton 
transverse momenta in the Breit frame, X^i PriPli the virtuality of the incident 
photon. In Fig. § the dependence of the two-jet cross section, in the cone scheme, on the 
renormalization and factorization scale factors S,r and E,f is shown. For scales related to 

pf (0 th^y are defined via 

= U (E PT{i)f , f^l = iF (E PT^)f ■ (3) 

i i 

For the two-jet inclusive cross section of Fig. the LO variation by a factor 1.43 is reduced 
to a 10% variation at NLO when both scales are varied simultaneously over the plotted 
range (solid curves). However, neither the LO nor the NLO curves show an extremum. The 
uncertainty from the variation of both scales for the NLO two-jet exclusive cross section in 
Fig. 2b (solid curves) is reduced to 5%. Furthermore, the two-jet exclusive cross section now 
has a maximum and is equal to the LO cross section for ^ = 0.5. Also shown is the ^ = C,r 
dependence of LO and NLO cross sections at fixed = 1/4 (dashed curves). In this case 
a maximum appears in the NLO inclusive and exclusive cross sections. However, the scale 
variation is stronger than in the ^ = = case. 

An alternative scale choice might he fij^ = fip = C, . The resulting ^ dependence is 
shown as the dotted lines for both the LO and NLO calculations. At LO the two scale choices 
give qualitatively similar results. However, with fij^ = fip = ^ Q^, the scale dependence does 
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not markedly improve at NLO. In addition a sizable ii'-factor is found, with K > 1 for small 
values of and K < 1 for very hard incident photons. We therefore use the jet transverse 
momenta in the Breit frame to set the scale and fix ^/j = = 1/4 in Eq. (^. A further 



discussion of the scale dependence of dijet cross sections at HERA can be found in Ref. [|13 

The effective i^-factors close to unity which are found in the cone and fc^ schemes could, 
in principle, be a coincidence arising from compensating effects in different phase space 
regions. It is important, therefore, to also compare LO and NLO distributions, in particular 
for those variables which define the acceptance region. 
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FIG. 3. Transverse momentum distribution in the lab frame for the jet with (a) minimal and 
(b) maximal transverse momentum. Results are shown for the two-jet inclusive cross section in the 
cone scheme in leading (dashed curves) and next-to- leading order (histograms). 



The transverse momentum distributions of the softest and the hardest jet in the labora- 
tory frame are shown in Fig. 3, for the cone scheme. The fairly small NLO corrections allow 
for reliable theoretical predictions. In general the largest radiative corrections are observed 
at small jet pt, as evidenced by the shape change in the PT,max distribution of Fig. 3b. 
The predictions are therefore expected to become more reliable for higher jet transverse 
momenta. A potential problem is the very steep p^^mm distribution, which, via the cut at 
5 GeV, introduces a strong sensitivity to the correct matching of the parton pt and the 
measured jet px- However, this is a general problem for all jet algorithms, i.e. the jet rate 
falls very rapidly as the required energy scale of the jets is increased. 

A more critical case is shown in Fig. 4 where, for jets defined in the /cr-scheme, the jet 
rapidity and the electron transverse momentum in the lab frame are shown. At NLO jets 
are produced somewhat more forward (in the proton direction) than at LO, see Fig. 4a. 
Hence, the rapidity cut at |?7j| = 3.5 has a stronger effect in NLO, which partially explains 
the relatively low i^-factor of 0.97 in the /c^-scheme. 



7 




FIG. 4. Rapidity distribution of the most forward jet (a) and transverse momentum distribu- 
tion of the scattered electron (b) in the lab frame. Results are shown for the scheme in leading 
(dashed curves) and next-to-leading order (histograms) for the two-jet inclusive cross section. 



Another observable which exhibits rather large NLO corrections is the electron transverse 
momentum distribution in Fig. 4b. The electron px becomes considerably softer at NLO, 
with an effective i^'-factor above unity at small Pt(^) and < 1 in the high transverse 
momentum region. In view of these shape changes the overall small change at NLO has to 
be considered a coincidence, tied to the choice of Pt(^) range. Since the electron transverse 
momentum and the of the event are very closely related, a similar change in the size of 
radiative corrections is obtained by choosing different bins. Very similar effects on the 
VmaxU) and Pt(0 distributions are also observed in the other jet definition schemes. As a 
result, a judicious choice of phase space region could generate very large or small i^-f actors 
which would indicate that, in these phase space regions, even the NLO calculation is fraught 
with large uncertainties. To avoid such potential problems, one should investigate the effect 
of the higher order corrections on those variables which are used to define kinematical cuts. 

A somewhat disturbing result of our calculation is the surprisingly large effective i^-factor 
in the W^-scheme, in particular since we disagree here with previous calculations The 
DISJET program 0, for example, gives a i^-factor very close to unity for a phase space 
region which is very similar^to the one considered above. The main difference to these earlier 



^About 10% are explained by differences in the virtual results for the "longitudinal projection" 
in the first paper of |^. This has been corrected in the latest version of 

^ The DISJET program does not allow to impose exactly the same kinematical cuts as discussed 
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calculations is the treatment of partons which are recombined into a jet. Previously, collinear 
and soft approximations were made for all partons which are coalesced into a single jet, i.e. 
for any pair of final state particles (including the remnant) with Sij < t/cut W^, neglecting 
terms of order t/cutW^. In our calculation, similar approximations are only made in a much 
smaller region, stj < Smin for any pair of initial and final state partons, and we are able to 
study the small Smin limit. Effectively, the approximations in Refs. correspond to the 
replacement of a massive two-parton system by a single massless jet. How well justified are 
these approximations? 

I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 400 I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 
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FIG. 5. Single jet mass effects at next-to-leading order, (a) Fraction of events in the cone 
scheme (solid curve), kT scheme (dashed curve), and Il^-scheme (histogram) with all jet mass 
to energy ratios below m/E, where E is the corresponding jet's energy in the parton center of 
mass frame. Negative values at small m/E are due to virtual contributions at m/E = 0. (b) 
Next-to- leading order transverse momentum (xt = Pt, solid histogram) and transverse-energy 
distribution (xt = Et, dotted curve) for the jet with largest pT and Et in the Breit frame, for the 
IF-scheme. The dashed curve shows the leading order result where both distributions are identical. 

In Fig. 5a the fraction of events is shown with at least one jet being more massive than 
m/E. Here m is the invariant mass and E the energy of the most massive of the jets in the 
parton center of mass frame. In LO m/E = since we are always using massless partons. At 



before for the W^-scheme; however, this is not relevant for the argument. 

^Our approximations in the soft and collinear region are stronger than the approximations made 
in [Qj5| where, effectively, also terms proportional to Sijlnsij are kept. These latter terms are 
neglected here. 
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NLO the median m/E is 0.44 in the VT-scheme {i.e. 50% of the NLO cross section hcs in the 
range m/E > 0.44), while substantially smaller values, 0.30 and 0.22 are found in the kx and 
cone schemes. The large single jet masses (compared to their energy) imply that coUinear or 
soft approximations for the partons inside a NLO jet are not generally permissable. Indeed 
the slow approach to the asymptotic region visible in Fig. 1 shows that only for two parton 
invariant masses squared Sij < 0.1 GeV^ are such approximations allowed. It is these mass 
effects which were not treated correctly in previous calculations. 

The very large median value of m/E in the ly-scheme implies that at NLO we are 
deahng with very different types of jets than at LO. At NLO at least one of the 1^-scheme 
jets extends over a large solid angle, it is a massive, slow moving object in the center of mass 
frame and, hence, very different from the pencil-like, massless objects called jets at LO. The 
typically small relativistic 7-factor of these jets has large kinematic effects. For example 
the difference between transverse energy and transverse momentum distributions of the jets, 
which are shown in Fig. 5b, becomes quite pronounced in the VF-scheme, an effect which is 
much smaller in the kx and cone schemes. 

The fact that very different kinematical regions are populated by the LO and NLO jets 
implies that radiative corrections may be very large as is indeed indicated by the large 
effective i^-f actor in the H^-scheme. The large corrections in turn imply that the predicted 
jet cross sections have large theoretical uncertainties and should not be used for precise QCD 
tests. Another measure of these uncertainties is the recombination scheme dependence of the 
predicted two-jet inclusive cross section. In the EO and P-schemes, where two partons are 
recombined to a massless jet, K is reduced to 1.41 and 1.29, respectively. This dependence 
first appears in the NLO calculation, where a jet may be composed of two partons. This 
internal jet structure, however, is only simulated at tree level and thus the dependence of 
the cross section on the recombination scheme is subject to potentially large higher order 
corrections. 

These uncertainties are small for jet algorithms with small recombination scheme depen- 
dence. In the JADE-algorithm the X-factor is reduced from 1.48 in the £'-schcmc to 1.36 
and 1.24 in the EO and P-schemes. For the cone (kx) scheme this recombination scheme 
dependence is reduced to the 3% (10%) level. The X-factor close to unity as well as this 
small recombination scheme dependence indicate that theoretical uncertainties due to higher 
order corrections are fairly small in the cone and kr scheme, and both schemes appear to 
be well suited for QCD studies at HERA. 
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